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Abstract 


Guerra, S. J., Lizárraga-Mendiola, L., £ Návar, J. (July- 
August, 2016). Riparian Communities along Longitudinal 
Gradients in Mexico's Northeastern San Juan River. Water 
Technology and Sciences (in Spanish), 7(4), 55-73. 


This research was conducted in three major tributaries 
of Mexico's northeast San Juan River with the major 
objectives of: (a) describing the diversity-abundance of 
riparian trees, benthic insects and fish faunal communities 
and (b) associating the fish and benthic insect communities 
to riparian tree communities, flow quality, and discharge 
parameters along a longitudinal gradient of water stress. 
Regardless of the high spatial variability, two gradients 
could be identified using multivariate analysis. Diversity of 
riparian trees, benthic insects and fish fauna was larger in 
more pristine and less perturbed upper in contrast to lower 
sampling sites. Bentic insect diversity was reduced and 
xerophitic riparian trees are colonizing lower sampling sites 
due in part to a larger biomass accumulation on rivers, higher 
stream temperatures and deteriorated river flow quality and 
diminished streamflow quantity. Benthic insects appear to be 
better adapted than fish fauna or riparian tree communities 
to river flow quality and quantity variability. The description 
and collection of riparian fish, benthic insects and riparian 
tree communities; as well as the ecological roles between 
biological groups are an important benchmark source of 
information and make this report useful to understand 
potential inherent sources of spatial and temporal variations 
of riparian communities. 


Keywords: Characterization of riparian habitats, alpha 
and beta diversity indices, longitudinal river segments, 
Multivariate analysis. 


Resumen 


Guerra, S. J., Lizárraga-Mendiola, L., € Návar, J. (julio-agosto, 
2016). Comunidades ripáricas en gradientes longitudinales en el río 
San Juan del noreste de México. Tecnología y Ciencias del Agua, 
7(4), 55-73. 


Esta investigación se realizó en tres principales tributaries del Río 
San Juan en el Nordeste de México, con los objetivos siguientes: (a) 
describir la diversidad-abundancia de las comunidades de árboles 
ripáricos, insectos bentónicos y peces y (b) asociar las comunidades 
de peces e insectos bentónicos a las comunidades de árboles ripáricos, 
los parámetros de la descarga y la calidad del agua de los ríos en 
gradientes longitudinales de perturbación del agua. Con el uso de 
análisis estadístico multivariados y a pesar de la alta variabilidad 
espacial se indentificaron dos gradientes principales. La diversidad 
de árboles ripáricos, insectos bentónicos y peces fue mayor en sitios 
más prístinos y menos perturbados en contraste con la menor 
diversidad encontrada en sitios aguas abajo con mayores índicios de 
alta perturbación. La diversidad de insectos bentónicos se redujo y 
árboles ripáricos xerofíticos se encuentran colonizando los sitios más 
perturbados, debido en partea la alta acumulación de materia orgánica 
en el agua, mayores temperaturas, calidad del agua deteriorada y 
menor descarga de ríos. Los insectos bentónicos como comunidades 
pioneras parecen estar mejor adaptados que los peces o los árboles 
riparios a las variaciones en la descarga y calidad del agua de los ríos. 
La descripción y la colección de los peces, insectos bentónicos y de 
árboles riparicos; asi como sus interacciones ecológicas entre grupos 
biológicos son una importante pieza de información y hacen que este 
reporte se útil para entender las fuentes inherentes potenciales de 
variación espacial y temporal de las comunidades ripáricas. 


Palabras clave: caracterización de hábitats ripáricos; índices de 
diversidad alfa y beta; segmentos longitudinales de ríos; análisis 
multivariados. 
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Introduction 


Riparian communities (e.g., fish fauna, benthic 
insects and riparian trees) typically contain 
diverse species assemblages that play a vital 
role in the health and productivity of aquatic 
ecosystems. Since rivers of semi-arid and sub- 
tropical landscapes are classified as intermit- 
tent and some segments as permanent streams 
(Levick et al., 2008), riparian communities adapt 
well to variable environments by developing 
resistance traits. High spatial and temporal river 
flow variability as well as a strong longitudi- 
nal zonation in temperature, discharge, water 
quality, substratum and biota feature semi-arid 
and sub-tropical rivers; with high flows quite 
often synchronized with seasonal precipitation 
events. In the dry season, low-velocity baseflow 
travel for short distances along the river chan- 
nel building waterholes and stagnant ponds 
that are often the only source of freshwater for 
fragmented riparian communities (Williams, 
2006). In the wet season, high flow or flood 
pulses re-connect river streams allowing sev- 
eral riparian faunal species to migrate locally 
and regionally to re-colonize neighboring and 
distant streams (Bunn, Davies, éz Winning, 2003; 
Nadeau éz Rains, 2007). The high spatial and 
temporal river flow variability, water quality, 
and predator-prey interactions make riparian 
species and communities to adapt well across 
these longitudinal gradients of river flows. 
However, resistance traits to river flow vari- 
ability must be different between biological 
groups (e.g., fish fauna and benthic insects) and 
linked to e.g., migration distances, river flow 
requirements, among others. 

Fish are ecologically important in rivers 
because they are often the apex predator in 
aquatic ecosystems. However, the diversity and 
abundance of fish species vary due in part to 
the patchy distribution of phreatic vegetation 
(Marsh-Mathews éz Mathews, 2000; Pusey éz 
Arthington, 2003) as well as due to the high 
temporal flow variability in semi-arid and sub- 
tropical rivers (Levick et al., 2008). Riparian 
vegetation is frequently correlated to fish faunal 


species assemblages as it provides (a) energy 
for primary and secondary productivity food 
webs for several terrestrial and aquatic biota 
(Pusey éz Arthington, 2003; Pettit et al., 2012); 
(b) food, cover, and refugee sites for aquatic as 
well as for terrestrial faunal species (Arthington, 
Balcombe, Wilson, Thomas, €z Marshall, 2005; 
Entrekin, Tank, Rosi-Marshall, Hollein, €z Lam- 
berti, 2009); (c) shade for moderating air and 
water temperatures (Bunn et al., 2003); (c) bank 
stability (Levick et al., 2008); and (d) filters for 
reducing contaminants, sediments, and debris 
(Roy, Rosemond, Leigh, Paul, 6: Wallace, 2003). 
Temporal stream flow variability, discon- 
tinuous riparian vegetation patches and fish 
faunal assemblages have also been linked to the 
diversity-abundance of benthic insects (Johnson 
£ Ringler, 2014); as these insects control nutrient 
cycles, primary productivity, decomposition and 
translocation of river materials (Peeters, Gylstra, 
éz Vos, 2004). Several benthic insect species are 
important consumers and important prey organ- 
isms (Whatley et al., 2014); controling the relative 
abundance of other faunal species in streams 
(Whatley et al., 2014). In summary, amongst the 
most important process shaped by the benthic 
insect community are: detrivory, grazing, and 
predator-prey interactions. Detrivory is the pro- 
cess of feeding on decaying organic matter and 
itis a major insect-mediated process. Benthic 
insects graze in nearly all aquatic ecosystems 
contributing largely to primary productivity 
of riparian ecosystems. Many benthic insects 
of all major groups function as predators but 
predatory fish often control the benthic insect 
community structure (Verdonschot, Didderen, 
€ Verdonschot, 2012; Whatley et al., 2014). 
Regardless of the importance of riparian 
ecosystems of semi-arid landscapes (Harper, 
Sanderson, €: McArthur, 1992; Levick ef al., 2008) 
some river components have not been properly 
described in Mexico's northeastern San Juan 
River (Návar, 2004, 2011). Important efforts had 
been previously carried out to describe the fish 
faunal species diversity along major tributaries 
of the San Juan River (Contreras-Balderas, 2000). 
However, more systematic studies are required 
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in order to understand sources of fish faunal as- 
semblage spatial variability across longitudinal 
gradients. There is also a lack of local informa- 
tion on the benthic insect community and how 
diversity and abundance of riparian fish fauna 
and benthic insect communities transition along 
longitudinal gradients of river flow as well as 
across different kinds of riparian vegetation. 
This information is critical to understand local 
controls of species diversity and abundance in 
order to eventually set appropriate restoration 
measures. Therefore, this report aimed: a) to 
describe the fish, benthic insect, and tree ripar- 
ian communities and b) to associate diversity 
and abundance of fish fauna and benthic insect 
communities to riparian trees, stream flow 
variability, river flow quality parameters and 
constituents for 10 sampling sites placed along 
three tributaries of Mexico's northeastern San 
Juan River. 


Materials and methods 
Description of the study area 


This research was conducted in Ramos and 
Pilon rivers, and the main stem of the San Juan 
River of Mexico's northeastern State of Nuevo 
Leon. These rivers rise in the upland watersheds 
of the Sierra Madre Oriental Mountain Range 
of Nuevo Leon and Coahuila. The main stem 
of the San Juan River flows through the Plains 
of the Northern Gulf of Mexico and drains into 
the lower Rio Bravo/Rio Grande at Camargo, 
Tamps., Mexico. "Ien sites were systematically 
chosen for sampling. When selecting sampling 
sites, the following sources of variation were 
considered: presence of riparian vegetation, 
water in the form of stream flow, waterholes 
or stagnant ponds; and a well defined channel. 
Sampling sites and major features of the land- 
scape are depicted in figure 1. 


Sampling 


Sampling consisted of measuring the physi- 
cal, chemical, and bacteriological parameters 


and constituents of streams, trapping fish and 
benthic insects, and measuring the dimensional 
characteristics of riparian trees. Sampling of 
fish, benthic insects, and riparian trees was car- 
ried out during the day from March to June of 
1999. At least one biological group was sampled 
weekly during the study period. Trees and 
shrubs were measured in two plots placed at 
the first terrace, immediately outside the river 
bank with dimensions of 40 m x 20 m. In each 
plot, we recorded diameter at breast height, 
dbh, top height, H, and canopy cover, Cc, for 
each individual tree. 

Nets dragged perpendicular to river flow for 
one hour collected fishes, taking care to drag on 
refugee sites beneath the river bank. Net was 5 
m in length and 1 m width. Surber-type sam- 
plers with 1 m in width collected benthic insects 
by dragging the bottom part of the stream across 
the channel or perpendicular to the river flow. 
Dragging was carried out three times at different 
places immediately below where the fish was 
collected. Collected macroinvertebrates were 
sorted and placed in large buckets filled with 
water and taken to the laboratory for counting 
and taxonomical identification. Collected fish 
and benthic insects were taxonomically identi- 
fied in Mexico's Laboratory of the Facultad de 
Ciencias Biologicas, UANL, in Monterrey, N.L., 
Mexico. 

Water samples were collected at midstream 
by way of wading and submerging plastic 
containers with a volume of 4 Lto a depth of 20 
to 30 cm, following the methodology proposed 
by Flores-Laureano and Návar (2002). Chemical 
constituents, and physical and bacteriological 
parameters were determined using laboratory 
methods of the Mexican Water Commission 
(CNA, 1993, 1997) following the standards of 
the Ministry of Urban Development and Ecol- 
ogy, denominated the NMX-AA-XX (Method 
Number)-YY (year) protocols for water analysis. 
For example, the method used for alkalinity de- 
termination is NMX-AA-36-1980. Measurements 
of stream velocity and river cross sectional area 
estimated discharge (Q) during each sampling 
event. Discharge measurements followed the 
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Figure 1. The location of the sampling sites along the San Juan River watershed in northeastern, Mexico. ChP = Cortijo- El 
Chapotal; BS] = La Boca; ChSJ] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda-La Corona; PSJ = El Porvenir; RR 
= Raíces; ES] = La Fragua; MS] = San Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; S] = San Juan. Note: Upstream sites 
are ChP, TR, RR, FR, HCP; downstream sampling sites are: BSJ, ChSJ, PSJ, FS], MS]. HCP was placed just before ites confluence 
with the major stem of the San Juan River. 


Mexican Water Commission standard as well. 
In addition, we recorded channel slope and lim- 
nological measurements for each sampling site. 


Data analysis 


We evaluated the alpha diversity by quantifying 
the species richness, Shannon, and Shannon's 
equitability parameters (Magurran, 2004). Shan- 
non's index accounts for both species abundance 
and species evenness and can be calculated by 
dividing H by H,.. (here H,.. =1n(S)). Equita- 
bility assumes a value between 0 and 1 with 
1 being complete evenness. Diversity indices 
were estimated for each biological group for 
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each sampling site. Equations described below 
evaluated each diversity index. The diversity 
index of Shannon - Weiner is: 


NT) 0 


Where p, = proportion of the total sample 
belonging to the ith species (1, /N), In = natural 
logarithm. The Shannon's equitability (E,,) is 
(Pielou, 1975): 


Es (Sp, “In(pi))/In(s) (2) 


Four diversity-abundance models fitted 
each sampling site for each biological group to 
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explain partition of limiting resources within the 
species and to draw conclusions on potential 
succession processes that shapes the commu- 
nity structure of each biological group. The 
geometric series, the log series, the truncated 
log normal distribution, and the broken stick 
diversity-abundance equations fitted the com- 
munity assemblage structure. The geometric 
series is: 


ni=N*Ck*k(1-k)" (3) 


Where ni = abundance of the ith species; 
N = total number of individuals or abundance; 
Ck = a constant, which assures that Eni = N; 
¡ = species rank. The logarithmic series is: 


ox,ax?/2,0x7/3,...,ax" /n (4) 


Where ax = the number of species with 
one individual; ax?/2 = the number of spe- 
cies with two individuals; a = (N(1-x))/x, and 
S/N = (1-x)/x(/In(1-x)). The truncated log 
normal series is: 


S(R) =S, exp(-aR)' (5) 


Where a = (20?) = a constant describ- 
ing the amount of spread in the distribu- 
tion; S(R) = number of species to occur in 
the Rth octave (=class) to the right or left of 
the symmetric curve; S, = number of spe- 
cies in the modal octave (the largest class). 
The broken stick model is (McArthur, 1988): 


S(n)=(S(S/1)/N)(1-1/N)* (6) 


Where S(n) = number of species in the 
abundance class that present n individuals; 
S = number of species; N = number of individuals 
or abundance. 

The similitude index of Morisita-Horn evalu- 
ated the beta diversity index, which is based on 
the richness and abundance of the most com- 
mon species. 

The index of Morisita-Horn, as described by 
Magurran (2004) is model (7): 


Ez 2) (an,bn,) 7) 
2 (da+dbjaNbN 


Where aN = total abundance of individuals 
in site A; LN = total abundance of individuals in 
site B; an, = abundance of the ¡-th species in site 
A; bn, = abundance of the ¡th species in site B: 


04 

an 

o 
Ñ aN? (8) 

bn? 
rs s 9 
AN? (9) 


The equal diversity null hypothesis between 
sampling sites was tested by the Shannon 
Weiner diversity index, since Hutcheson (1970) 
described the t test given by equation (10): 


H,-B, 


f= 
(VarH, +VarH,)” 


(10) 


Where H, ,= indices of Shannon in samples 
l and 2; Var H, and H,= deviance of the index 
of Shannon in samples 1 and 2 with degrees of 
freedom given by equation (11): 


(VarH, +VarH As 


(1) 
(VarH, ) /N, +(VarH,)' /N, 


df = 


Where N, and N, = number of individuals in 
samples 1 and 2. 

Detrended canonical correspondence analy- 
sis was performed for each biological group us- 
ing river flow discharge and quality parameters 
as the environmental variables to understand 
clusters and patterns of species and sampling 
sites as well as the major sources of variation 
explaining clusters and patterns. In addition, 
several regression equations were tested to 
predict e.g., fish fauna species assemblages as a 
function of the riparian tree or the benthic insect 
species as exogenous variables using several 
alpha diversity indices. 
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Results total dissolved solids, and sulfates (table 1). 

Tributaries El Pilon (ChP, HCP) and Ramos (FR, 
Stream water parameters and constituents RR, TR) had the fewest number of exccedances 

in contrast to the main stem of the San Juan Riv- 
Several constituent concentrations and param- er (BSJ, CHSJ, PSJ, FSJ, and MSJ). Water quality 
eters exceeded the US-EPA National Primary parameters and constituents deteriorated from 
Standards for drinking water for chloride, fecal upper (ChP and FR) to middle stream (ChSJ, 
coliform bacteria, total coliform bacteria, color, HCP, and PSJ) and to lower sites (FS] and MSJ). 


Table 1. Concentration statistics for measured constituents and parameters of stream water quality of the San Juan River, N.L., 


Mexico. 
Site 
Parameters 

ChP BSJ ChSJ FR HCP PSJ RR FS] MSJ TR 

Total acidity (mg 1* CaCo,) 0 6 7 2 5 0 0 2 0 0 

Phenolic alkalinity 7 4.5 0 0 0 3 6 0 1 1 

Total alkalinity () 150 228 275 252 130 157 167 366 141 

Chloride (mgl*) 1 52 61 43 66 E 4 79 17 

Fecal coliform B. (NMP 100 ml") 20 5M 24M 8 000 2 20 2 170 700 

Total coliform B. (NMP 100 ml”) 90 5M 24M 8 000 60 20 |90000 | 17000 700 

Color, Upt-Co 5 41 36 46 152 5 8 36 12 

Electric conductivity (umos cm”) 601 931 1095 1558 956 357 646 1065 378 

DBO (mgl”) 2 5.6 3 4 2 2 2 2 2 

DQO (mgl”) 3 25 19 18 9 5 5 6 5 

Total hardness (mgl* CaCO,) 324 364 317 690 329 184 277 447 182 

Ca hardness 261 301 296 518 165 156 211 156 85 

Mg hardness 63 63 21 172 164 28 66 291 97 

Oil content 0.47 6 2 2 1.31 0.69 0.91 6.03 0.7 

Nitrite-nitrate 0.002 | 0.011 0.012 0.016 0.003 | 0.001 0.004 0.034 | 0.003 

pH 7.88 8.02 7.84 7.78 7.82 UI 7.46 7.84 7.78 

e SAAM (mgl”) 0.02 4.3 0.021 0.035 0.422 | 0.011 0.037 0.07 | 0.058 

e T. solids (mgl”) 473 692 806 1319 753 287 459 808 265 

S T. suspended solids 3 29.2 38 61 66 2 4 11 3 

3 T. disolved solids 470 | 664.3 768 1258 687 235 455 797 262 

E S. solids 0.1 4.3 0.1 0.1 0.3 0.1 0.1 0.1 0.1 

So Sulfate (mg]l”) 148 189 251 590 266 36 115 197 36 

Z SW temperature (C) 23 23 23.5 28.5 31 25 29 30 34 

Y Env. temperature ("C) 24.5 25 24.5 28.5 34 33 30 27 30 

E Turbidity (NTU) 0 5 7 0 7 28 0 1 9 2 

=] Channel depth(m) 0.37 | 0.34 0.75 0.32 1 0.5 0.33 0.55 0.90 0.80 

? Channel width(m) 12.6 9.1 9.0 11.2 5 5 16.4 7 20 12.5 

= Channel area(m?) 4.7 3.0 6.8 3.6 5.0 2.5 5.4 3.9 18.0 10.0 

S Discharge (més”) 0.22 1.20 1.30 0.90 1.20 0.23 0.07 0.30 0.24 0.11 
E Bold and Italic values of constituents and parameters exceed the US-EPA National Primary Standards for drinking water. ChP = Cortijo, El Chapotal; 
= BSJ = La Boca; ChS] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; RR = Raíces; FS] = La Fragua; 
Y MS] = San Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; S] = San Juan. T = total; Env. environmental; SW = surface water; S = suspended. 
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A principal component analysis reduced 
the multiple dimensions of the water quality 
parameters and constituents and indicated that 
92% of the total variance was explained by the 
first seven components. Therefore variables 
with the largest weight on each component 
were selected to explain the intrinsic sources of 
variability. These included phenolic alkalinity, 
electric conductivity, oil content, nitrite-nitrate 
concentration, total and dissolved solids, stream 
discharge, temperature, and turbidity. 


The riparian tree community 


Tree communities of Mexico's northeastern San 
Juan River are diverse and abundant (table 2). 
We recorded a total of 788 trees of 24 genera, 
16 families, and 13 orders. Average stand den- 
sity was only 500 trees ha”, with a standard 
deviation of 488 trees ha” per sampling stand. 
Taxodium mucronatum, atypical riparian tree and 
an abundant tree species, was absent in four out 
of ten downstream sites (PSJ, RR, FSJ, and MSJ). 


Table 2. Diversity-abundance of riparian tree communities in 10 sites placed along the San Juan River of northeastern Mexico. 


Species pS Abundance 
ChP | BS] | ChSJ]| FR | HCP | PS] | RR | FS] | MS] | TR 

Taxodium mucronatum Ten. 91 31 27 76 12 4 241 
Platanus occidentalis L. 66 6 31 29 1 133 
Salix nigra Marsh 9 4 2 1 12 6 10 129 
Baccharis glutinosa (R £P) Pers. 72 14 7 93 
Acacia farnesiana (L.) Willd 6 2 2 1 5 9 19 44 
Parkinsonia aculeata L. 18 18 
Ig (Terán « Berl.) 2 16 de 
nas ebano (Berlandier) C.H. 1 12 3 1 7 
Cephalanthus occidentalis L. 1 1 4 1 2 7 16 
Melia azederach L. 1 3 6 3 13 
Leucaena leucocephala Benth 10 10 
Havardia pallens Standl. 2 2 1 1 3 9 
Dodonea viscosa Jacq. 9 
Morus microphylla L. Y 1 8 
Fraxinus sp. Gray. 3 1 2 6 
Juglans sp. Engelin. 1 1 1 2 5 
Celtis laevigata Willd 1 1 2 1 5 
Celtis pallida Torr 2 2 4 
Populus tremuloides Mich. 3 3 
Caesalpinia mexicana Gray 1 1 2 
Zanthoxylum fagara (L.) Sarg 2 2 
Prosopis laevigata (Willd.) H.C. 

Johnst. E E 
Sargentia greggii Wats 1 1 
Diospyros texana Scheele. 1 1 
Abundance 270 60 72 142 18 90 20 62 20 34 788 
Richness 10 9 7 14 5 4 7 8 4 7 


ChP = Cortijo, El Chapotal; BS] = La Boca; ChS] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; RR 
= Raíces; FS] = La Fragua; MS] = San Miguelito; TR = Las Trancas; P = Pilon; R = Ramos; S] = San Juan. 
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T. mucronatum was also quite rare along sam- occidentalis, P. tremuloides, S. nigra) in most up- 
pling site TR with only four trees in two plots per and middle sampling sites for fitting any 
of 800 m? each. Riparin tree communities with diversity-abundance model (table 3). 
this typical phreatophitic tree as the dominant The beta diversity described by the Morisita- 
tree species are located in upper stands of main Horn index was quite small for all sampling 
tributaries. A mixture of riparian and xerophytic sites, with an average (confidence interval) of 
tree species dominates lower riparian corridors. 0.28 (+ 0.09). That is, the riparian tree commu- 
In lower sampling stands, we measured a nity is quite different between sampling sites 
large number of xerophitic tree species. Acacia with only a few species colonizing neighboring 
farnesiana was frequently observed in 7 out of stands. In general, sampling sites placed at a 
10 sampling sands, although it accounted for similar relative location, e.g., the upper sites, 
less than 10% of the total abundance. Other recorded large beta indices (ChP, BSJ, ChSJ, FR). 
xerophitic tree species found preferentially in The beta diversity dropped sharply at lower 
lower stands (P. aculeata, E. anacua, E. ebano, L. sites of the same tributary. 
leucocephala, H. pallens, C. laevigata, C. pallida, The large variability of the riparian tree as- 
Z. fagara, P. laevigata, S. greggíi, and D. texana) semblages makes to fit well all successional dis- 
accounted for less than 50% of the recorded tributions: the logarithmic series (intermediate); 
diversity. the log-normal (intermediate); the geometric 
Tree diversity varied between sampling (pioneer) or the broken stick (late successional) 
sites. The t test on the Shannon-Weiner diversity models. 
index showed tree diversity was smaller at the The detrended canonical correspondence 
middle stream, PSJ, sampling site since it had analysis (DCCA) indicated the first two tree 
only 4 tree species and 90 trees ha”. Evenness species axis accounted for 35% of the total 
is larger at lower (HCP, FSJ, MSJ, and TR) than species deviance and 60% of the total tree 
at the upper sampling sites (ChP, BSJ, and FR). species-environment relation. Two longitudinal 
All diversity-abundance models predicted well gradients could be noted in figure 2. 
the community structure pointing at a succes- The first one moves along the Y axis from 
sional equilibrium of the plant community. the center (ChSJ, ChP), to the middle (FR, HCP), 
In spite of the presence of xerophitic species, to the upper (RR) part of the Y axis. Tree spe- 
the tree riparian community has sufficient late cies Taxodium mucronatum, Populus tremuloides, 
successional tree species (T. mucronatum, P. Juglans spp., Morus microphylla, and Leucaena 


Table 3. Diversity indices and diversity-abundance models fitted to the riparian tree community of the San Juan River of 


ChP = Cortijo, El Chapotal; BS] = La Boca; ChS] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; RR = 
Raíces; FS] = La Fragua; MS] = San Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; S] = San Juan. HO and HA = null and alternate hypotheses. 
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Figure 2. Results of the DCCA analysis on the riparian community of trees of the San Juan River, N.L., Mexico. ChP = Cortijo, El 
Chapotal; BS] = La Boca; ChS] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; RR = 
Raíces; FS] = La Fragua; MS] = San Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; S] = San Juan. 


leucocephalla clustered around upper sampling 
sites (FR, ChP, ChSJ, RR, BSJ) featured by low 
stream temperature, high discharge, best water 
quality parameters and constituents, and small 
channel width, area and depth (figure 2). The 
second longitudinal gradient is located along 
the X axis, moving from the center (BSJ, ChSJ, 
ChP) to the right (MSJ, FSJ, PSJ) of figure 2. 
The tree species, Platanus occidentalis, Dodonea 
viscosa, Salix nigra, Bacharis spp., and Caesalpinia 
mexicana clustered around lower sites (ChSJ, FSJ, 
PSJ, and MSJ). Lower sampling sites (FSJ, PSJ, 
MSJ, and TR) are characterized by increased 


stream flow temperature and pollution levels, 
channel depth, and reduced discharge. Acacia 
farnesiana, Ehretia anacua, Diospyros texana, Pithe- 
cellobium pallens, and Fraxinus spp. agglomer- 
ated around middle and lower sites (TR, MSJ, 
FSJ, PSJ). 


The fish community 


The San Juan River hosts an important fish 
diversity and abundance (table 4), since we 
captured a total of 5 608 fishes classified in 17 
genera and 8 families. Gambusia affinis, a grazer 
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Table 4. Diversity-abundance of fish species captured in 10 sites of three tributaries of the San Juan River, Nuevo Leon, Mexico. 


Code Species E 
ChP | BSJ FR PSJ RR FS] TR | Abundance 
S  |Gambusia affinis 2 8 1 2 260 | 1892 2580 
S  |Poecilia formosa 8 227] 58 1 414 DO 809 
S Astyanax mexicanus 4 2 1 392 130 635 
S  |Dionda couchi 240 14 176 430 
AS  |Notropis amabilis 224 128 22 30 7 411 
S  |[Campostoma anomalum 50 132 187 369 
S  |[Poecilia mexicana 16 66 20 102 
E [Tilapia spp. 48 12 60 
S  |Micropterus salmoides 1 38 4 11 54 
S  |Cichlasoma (Herichtys) cyanoguttatum 2 2 a 1 20 8 42 
S  |Cyprinella lutrense 37 37 
S  |Lepomis megalotis 8 18 2 4 32 
S  |Cyprinella rutila 17 17 
AS  |Moxostoma congestum 15 1 16 
V  |Etheostoma grahami 8 8 
S  |Lepomis macrochirus 3 3 
V  |lctalurus lupus 1 1 2 
NOM |Heterandria cf. Bimaculata 1 1 
Abundance 270 284 528 148 109 298 | 1366 | 2605 5 608 
Richness 5 7 5 3 7 12 9 12 


ChP = Cortijo, El Chapotal; BS] = La Boca; ChS] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; RR 
= Raíces; FS] = La Fragua; MS] = San Miguelito; TR = Las Trancas; P = Pilon; R = Ramos; SJ = San Juan. S = secure, AS = apparently secure, V = 


imperiled, NOM = NOM-059-SEMARNAT-2001. 


fish species that colonizes the lower parts of 
rivers accounted for close to 50% of the total fish 
abundance. Of the 18 fish species sampled, 7 are 
grazers, 3 are carivores, 5 are invertivores, and 
the remaining 3 species have unknown feeding 
habits. Gambusia affinis, Poecilia formosa, and 
Cichlasoma cyanoguttatum have large ecological 
plasticity or widespread habitat since these spe- 
cies were captured at 7 out of 8 sampling sites. 
On the other side, the fish species L. macrochirius, 
E. grahami, H. jonesi, and C. lutrense have a more 
restricted longitudinal habitat since they were 
captured at only one site. 

Most captured fish species are classified as 
secure in terms of the round global status (Na- 
ture Server Explorer, 2008). However, two fish 
species are classified as apparently secure (AS) 
in southern USA (Notropis amabilis and Moxos- 


toma congestum). Etheostoma grahami, Ictalurus 
lupus, and Compostoma anomalum are endemic 
species to Mexico's northeastern San Juan River 
and southern Texas. The first species is locally 
classified in the NOM-059-SEMARNAT-2001, 
and it was captured only at RR, a sampling site 
with the best river flow quality parameters. 
Heterandria cf. jonesi is an invasive species first 
recorded during this study in the San Juan 
River watershed, since it was absent on the list 
reported by Contreras-Balderas (2000) and it 
was captured only at ChSJ, a heavily polluted 
sampling site. 

The fish diversity indices and diversity- 
abundance models are reported in table 5. 
Sampling sites placed at lower sampling sites 
such as the Ramos River (RR and TR) recorded 
the largest fish diversity with 12 out of the 18 
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Table 5. Diversity indices and diversity-abundance models fitted to the fish communities of the San Juan River in northeastern 


Mexico. 
Diversity pe 

ChP BSJ ChSJ FR PS] RR FSJ TR 
Shannon 0.48 0.82 0.58 0.44 1.13 1.38 1.69 1.00 
Shannon equitabilty E,, 0.30 0.42 0.36 0.40 0.58 0.56 0.77 0.40 
Models 
Geometric series HA HA HA HA HA HA HA HA 
Logarithmic series HO HO HO HO HO HO HA HO 
Log-normal HO HO HO HO HO HO HO HO 
Broken stick HA HA HA HA HA HO HA HA 


ChP = Cortijo, El Chapotal; BS] = La Boca; ChSJ] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; RR = 
Raíces; FS] = La Fragua; MS] = San Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; S] = San Juan; HO and HA = null and alternate hypotheses. 


species. Sampling sites placed along the upper 
Ramos River (FR) and Pilon River (ChP) had the 
lowest fish diversity-abundance. 

The Shannon index varies notoriously be- 
tween sampling sites. The fish diversity was 
lowest at the upper sites (ChP, ChSJ, and FR) 
and largest at the lower sites FS] and RR. FSJ 
and RR recorded one of the largest abundance 
and species richness of all sampling sites. The 
rest of the sampling sites differ in a significant 
manner from these two groups. Evenness had 
a tendency to increase from upper to lower 
sampling sites. The log normal and logarithmic 
series diversity-abundance models fitted the 
fish community structure better. The pioneer 
and late succesional models represented by the 
geometric series and the broken stick models 
did not fit the fish assemblage structure for any 
particular site. Since the geometric series and 
the broken stick model apply well to pioneer 
and late successional stages or well balanced 
communities at any developmental stage, 
recorded data shows and fitted models predict 
the fish community is in the secondary stage of 
succession. 

The beta diversity index was also quite 
small for all sampling sites, with an average 
(confidence interval) of 0.18 (+ 0.12). That is, 
the fish community differs between sampling 
sites with only a few species distributing and 
colonizing neighboring sampling sites. Similar 


dominant fish species were observed at sites of 
similar location; e.g., upper sites ChP and RR 
and middle sites BS] and FR. Sampling sites 
ChSJ and TR recorded the largest beta index 
(> 0.85). However, fish fauna collected in the 
same tributary (e.2., FR, RR, TR; or ChP, HCP) 
had small beta diversity indicating a lack of 
connectivity between fish species. 

The fish assemblages fitted better the loga- 
rithmic series and the log-normal model unlike 
the geometric (pioneer) or the broken stick (late) 
models stressing the community is found within 
the intermediate stages of succession. 

The DCCA analysis carried out on the loga- 
rithmic transformation of the abundance data 
showed the first two species axis accounted for 
47% of the total fish species variance and 35% 
of the total fish species-environment relation 
(figure 3). 

Two longitudinal gradients were noted in 
the fish structure assemblages. The first gradient 
moves along the Y axis and is described by E. 
grahami, C. rutila, M. congestum, N. amabilis, D. 
melanops, N. jamezanus, and 1. lupus, as well as by 
sampling upper sampling sites ChP, RR, FR, BSJ. 
Low streamflow temperature, high discharge 
rate, as well as low nitrite-nitrate concentra- 
tions, suspended solids, and total suspended 
solids characterize these places (figure 3). The 
second longitudinal gradient moves parallel 
to the first one but on the right hand side of 
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Figure 3. Results of the DCCA analysis on the fish community of the San Juan River, N.L., Mexico (TD Solids = total dissolved 
solids, Feno Alkalinity = fenolic alkalinity, T. Solids = total solids). ChP = Cortijo- El Chapotal; BS] = La Boca; ChSJ = Ejido 
Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PSJ = El Porvenir; RR = Raíces; FS] = La Fragua; MS] ='San 
Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; SJ = San Juan. 


figure 3. The fish species C. lutrense, P. spp., G. 
affinix, A. mexicanus, Tilapia spp., P. formosa, C. 
cyanoguttatum, L. macrochirius, and M. salmonoids 
colonized the lotic waterholes and ponds of 
lower sites FSJ, PSJ, TR featured by increased 
streamflow temperature, small discharge rate, 
larger nitrite-nitrate concentrations and sus- 
pended sediments. 


The benthic insect community 


The benthic insect community was the most 
diverse biological group found on the Mexico's 
northeastern San Juan River (table 6). The 
data source was made of 865 individuals, which 
belong to 39 different species, 34 genera, 23 
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families, and 6 orders. Families Ephemeroptera 
and Coleoptera contained the highest species di- 
versity with 38% and 37%, of the total diversity, 
respectively. The family Trichoptera recorded the 
least relative abundance, with only 3.4% of the 
total diversity. Middle sampling sites showed 
the largest benthic insect diversity (PSJ) ac- 
counting for approximately 50% of the total 
observed diversity. However, the upper, most 
pristine sampling site along the Ramos River 
recorded the largest diversity of rare species. 
The Ramos River had the largest benthic insect 
diversity with 12 out of the 18 species observed 
here. Sampling sites distributed along the upper 
Ramos and Pilon rivers had the smallest benthic 
insect diversity and abundance. 


Guerra et al., Riparian Communities along Longitudinal Gradients in Mexico's Northeastern San Juan River 


Table 6. Diversity-abundance of benthic insect species captured in 10 sites of three tributaries of the San Juan River, 


Nuevo Leon, Mexico. 
Code Species Sa Abundance 
ChP | BSJ |[ChS]| FR | HCP | PS] | RR | FS] | MS] | TR 
Ps  |Psephenus sp. 104 1 2 9 87 86 289 
Ca  |Caenis sp. 27 | 104 131 
Tr  |[Tricorythodes sp. 1 1 27 10 37 3 79 
Te  |Tenagobia sp. se 55 
Ar  |Argía sp. 1 1 14 9 4 3 10 48 
Ba  |Baetis sp. 2 Y) 1 8 3 11 13 47 
Le  |Leptophlebia sp. 30 1 31 
Ch  [Choroterpes sp. 1 22 5 1 1 30 
An  |Ancyronyx sp. 5 16 21 
Ci |Chironomus sp. 2 22 il 17 
He  |Helicopsyche sp. 8 5 13 
Am  |Ambrysus sp. O 1 2 2 2 12 
Ce  |[Cernotina sp 12 12 
Cl [Callibaetis sp. 9 9 
Ta  [Tanytarsus sp. 9 
La  |Larsia sp. 1 1 4 1 7 
Sc  |Scirtidae 1 5 6 
Ma  |Macrothemis sp. 5 5 
Co  |Corixidae 1 1 2 4 
Tn  [Tanypodinae 4 4 
Ap  |Anopheles pseudopunctipennis 3 1 4 
Ph  |Phyllogomphoides sp. 3 3 
Ai  |[Anisoptera 1 2 3 
Ly  |Lype sp. 1 2 3 
Lb  |Labrundinia sp. 2 1 3 
Ne  |[Neurocordulia sp. 1 1 2 
Be  |Berosus sp. 1 1 2 o 
Br  |Brachycentrus sp. 22 2 ne 
Th  |[Thienemannimyia sp. 2 2 = 
Ab  |Ablasbemyia sp. 2 2 S 
Ce  [Chironomus? 2 2 
Is  |Isonychia sp. 1 1 E 
Lp  |Leptohyphes sp. 1 1 E 
Ht  |Hetaerina sp. 1 1 3 
Ag  |Argiallagma sp. 1 1 _ 
Hi  [Hidrophilidae 1 1 3 
Po  |Polypedilum sp. 1 1 5 
Ni [Nímbocera sp. 1 1 E 
Lu  |Leucotabanus sp. 1 1 E 
Abundance 14 8 | 123 | 140 16 33 23 | 171 | 212 [| 125 865 3 
Richness 7 2 7 | 18 5| 12 7 9| 13 9 3 
S 
ChP = Cortijo- El Chapotal; BS] = La Boca; ChSJ] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; RR Ed 
= Raíces; FS] = La Fragua; MS] = San Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; S] = San Juan. E 
E 
E 
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Most frequently observed benthic insect appears that resources are well distributed 
species (Argia spp. (8), Baetis spp. (8), Psephenus between the few species observed. 
spp. (6) and Tricorythodes spp. (6). Macrothemis The beta diversity index was quite small for 
spp., Callibaetis spp., Tanytarsus spp. and all sampling sites, with an average (confidence 
Phyllogomphoides spp.) are classified as rare interval) of 0.22 (+ 0.08). That is, the benthic 
species in these rivers of northeastern Mexico. insect community diversity varies greatly be- 
In this riparian ecosystem, diversity generally tween sampling sites. The probability of finding 
increased with abundance. The least polluted values > 0.75 was similar between rivers and 
Ramos River reported the highest insect diver- quite low in neighboring sampling sites and the 


probability diminished even more in sampling 
sites farther appart. That is, there is a lack of 
spatial tendencies or patterns of species con- 
nectivity in the beta diversity index between 
and among sampling sites. The patchy structure 
of benthic communities given by the small beta 
index shows the existence of place to place dif- 
ferences in environmental features such as food 
distribution, substrate, discharge, water quality, 
among others. Stream disturbances may also be 
controlling the lack of connectivity patterns in 
community diversity structure. 

The large variability of the benthic insect 
assemblages makes to fit well all successional 
distributions as well: the logarithmic series 
(intermediate); the log-normal (intermediate); 
the geometric (pioneer) or the broken stick (late 
successional) models. 


sity and one of the largest abundances with 
greater than 30% of the total number of benthic 
insects captured. 

The Shannon-Weiner diversity index for 
benthic insects varied between sampling 
sites. Sampling sites BS] and ChSJ recorded 
the smallest and sampling sites FR and PSJ 
recorded the largest alpha diversity index. 
Evenness showed a tendency to decrease from 
upper to lower sampling sites. The logarithmic 
series model fitted the benthic insect diversity- 
abundance structure for 9 out of 10 sampling 
sites. This community is in a secondary stage of 


succession. Sampling sites with intermmediate 
diversity-abundance parameters (FSJ, MSJ, TR) 
had a tendency to show species of the second- 
ary stages of succession. Sampling sites with 


the largest Shannon diversity index (FR and The first two DCCA axes accounted for 44% 
PSJ) as well as with the highest water quality of the total benthic insect species variance and 
parameters tended to be fit by all diversity- 45% of the total benthic insect species-environ- 
abundance models. That is, the abundance ment relation. Two longitudinal gradients could 
distribution of species is well balanced and be identified in figure 4. 


Table 7. Diversity indices and diversity-abundance models fitted to the benthic insect communities of the San Juan River 


ChP = Cortijo- El Chapotal; BS] = La Boca; ChS] = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; 
RR = Raíces; FS] = La Fragua; MS] = San Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; S] = San Juan. HO and HA = Null and alternate 
hypotheses. 
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A 

[7 

a 

Ss in northeastern Mexico. 
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E Site 

E Diversity 

ES ChP BSJ ChS] FR HCP PSJ RR FS] MS] TR 
Z Shannon 1.73 0.38 0.58 2.33 1.24 2.16 1.63 1.38 1.55 1.12 
+ Shannon equitabilty E,, 0.89 0.55 0.30 0.81 0.77 0.87 0.84 0.63 0.60 0.51 
3 Models 
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Figure 4. Results of the DCCA analysis on the benthic insect community of the San Juan River, N.L., Mexico (TD Solids = total 
dissolved solids, Feno Alkalinity = fenolic alkalinity, T. Solids = total solids). ChP = Cortijo-El Chapotal; BS] = La Boca; 
ChSJ = Ejido Chihuahua; FR = El Fraile-Terreros; HCP = Hacienda, La Corona; PS] = El Porvenir; RR = Raíces; FS] = La Fragua; 
MSJ = San Miguelito; TR = Las Trancas; P = Pilón; R = Ramos; S] = San Juan. 


The first one is characterized by the benthic 
insect families Ambrysus, Helicopsyche, Lype, 
Psephenus, Leucotabanus, Argia, and Ancyronyx 
colonizing sampling sites FR, ChP, RR, TR, and 
ChS]J. These families inhabit mostly upper tribu- 
taries featured by the highest stream water qual- 
ity indicators, with low concentrations of sus- 
pended solids, nitrite-nitrate, oil content, fenolic 
alkalinity; as well as with lowest streamflow 
temperatures and highest discharge rates. The 
second longitudinal gradient is characterized by 
the benthic insect families Neurocordulia, Berosus, 
Baetis, Arisoptera, Larsia, Caenis, Ablasbemyia, 
Nimbocera, Chironomus, Hidrophilidae, Anopheles, 
and Tenagobía distributed along sampling sites 
ChP, PSJ, BSJ, MSJ, and HCP. 


Discussion 


The studied riparian communities (trees, fish 
and benthic insects) could be consistently clas- 


sified by the DCCA analysis into upper and 
lower river-segment communities for all three 
major biological groups. A lack of connectiv- 
ity of these biological groups between upper 
and lower sampling sites was identified as the 
major source of variation by several indicators; 
the DCCA analysis, the Shannon-Wiener alpha 
diversity index, and the Morisita-Horn beta 
diversity index. 

Springs feeding on river channels in upland 
ridges along the piedmont characterize upper 
river segments; with cold flow temperatures, 
high water quality indicators, high flow, and 
a well defined channel cross section. Typical 
riparian tree communities (T. mucronatum, P. 
occidentales, Juglans spp., Populus spp.) colo- 
nize these stands. Fish fauna assemblages (E. 
grahami, C. rutila, M. congestum, N. amabilis, 
IL. lupus, and N. jamezanus) and benthic insect 
families (Ambrysus, Helicopsyche, Lype, Psephe- 
nus, Leucotabanus, Argía, and Ancyronyx) thrive 
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on these rivers and they are indicators of good 
water quality. 

Stream temperature and pollution levels 
increase and flow diminishes at lower river 
segments. Diversity of typical riparian tree 
communities diminishes from upper to lower 
sampling sites. This tendency is related to the 
dissapearance of riparian phreatophitic and 
the appearance of xerophitic tree species. Other 
sources of variation such as soil type, altitude 
above sea level, and geology may be controling 
the presence of T. mucronatum across longitu- 
dinal riparian gradients (Hudak é Ketcheson, 
1992; Rosentreter, 1992). The disappearance of 
typical riparian trees as the river flows away the 
piedmont can be one of the adjustments riparian 
vegetation use in response to the modification of 
river flow parameters (Johnson, 1998). 

The fish diversity has a slight tendency to 
decline from upper to lower river segments. G. 
affinis, Poecilia spp., C. lutrense, A. mexicanus, L. 
megalotis and Tilapia spp. colonize well lower 
and middle sites (TR, FSJ, and ChSJ). These spe- 
cies assemblages are characterized by increased 
pollution levels, flow temperature and reduced 
discharge. Several chemical constituents (deter- 
gents and sulfates) and bacteriologic parameters 
(fecal coliform bacteria) exceed the NOM-001- 
ECOL-1996 standard that sets the maximum 
permisible limits on water quality parameters 
as well as the standard for aquatic organisms 
CE-CCA-001/89. De Lange (1994) and Villarreal 
(1983) noted several fish species described im- 
mediately above are indicators of poor stream 
water quality with depleted oxygen, high sedi- 
mentation rates, and degraded canopy cover. 

The benthic insect communities thriving 
along lower river segments (Neurocordulia, 
Berosus, Baetis, Arisoptera, Larsia, Caenis, Ablasbe- 
myia, Nimbocera, Chironomus, Hidrophilidae, 
Anopheles, and Tenagobia) have been described 
as good indicators of deteriorated stream water 
quality (Phillips € Rainbow, 1993; Merrit éz 
Cummins, 1996), high stream temperature and 
low flows (Haidekker €: Hering, 2007). The pres- 
ence of these species in polluted streams relates 
to their adaptation to these conditions as well as 
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to the hability to degrade organic matter since 
Leyton and Voshell (1991) found these families 
can recycle up to 20% of the organic matter in 
streams. On the other side, Lenat (1993) noted 
the family Tricorytoides has the potential to be an 
indicator of good water quality. 

River fragmentation by one of the most acute 
drought spells of the last century (1990's) may 
partially explain the lack of hydrological con- 
nectivity between upper and lower fish fauna 
and benthic insect communities. The decadal 
drought spell that affected Mexico's northern 
watersheds recorded consistently below average 
precipitation for most of the decade and started 
sometimes during 1989-1990 and finished some- 
times during 2002-2003 (Návar, 2015). During 
this period, river discharge stopped flowing 
in numerous occasions in several gauging sta- 
tions along perennial tributaries of the San Juan 
River (Návar, 2011). In spite of this acute natural 
disturbance that restricted long-range migration 
of several fish species (Ictalurus lupus, Lepomis 
megalotis, and Moxostoma congestun) as well as of 
several benthic insect species (Leptophlebia spp.) 
in upper river segements, the trophic structural 
levels are in good condition as the ecosystem 
contains primary consumers (Helicopysche spp. 
consumes algae), predators (Notropis amabilis) 
and preys (Chironomus spp.). The trophic chains 
are also well defined in lower river, polluted 
segments by the presence of predators (Ictalurus 
lupus and Gambusia affinis), preys (Anopheles spp. 
and Chironomus spp.) and the largest predator, 
the exotic fish species Tilapia spp. 

Fish fauna and benthic insect communities 
are then well adapted to drought episodes of 
several time-scales by developing resistance 
traits. Adaptation traits appear are more vari- 
able between the benthic insect communities 
than between the fish fauna communities as the 
fish longitudinal gradient are better identified 
than the benthic insect longitudinal gradients. 
Fish is adapated to migrate short (Lepomis 
megalotis, Gambusia affinix, Poecilia spp.) and 
long-range (Ictalurus lupus, Ictalurus spp., Tilapia 
spp., Notropis amabilis, Moxostoma congestum) 
distances but they require of continuous river 
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flow connectivity to colonize preferentially up- 
per river segments. Benthic insects migrate but 
only for short distances and they may not need 
constant river flow connectivity all the time to 
meet the physiological needs. 

The breaking of two fish faunal longitudinal 
gradients may partially be also a response to 
other sources of disturbance such as the pres- 
ence of the introduced exotic fish Tilapia spp. 
This large predator was restricted predominant- 
ly to lower river segments characterized by high 
flow temperature and increased pollution levels. 
A lack of river connectivity during seasonal, 
inter-annual or decadal dry spells may modify 
local upper fish community assemblages as 
these conditions restrict migration of large fish 
predators such as Tilapia spp. and Ictalurus spp. 


Observations between communities 


The tree, fish, and benthic insect species diver- 
sity described by the E,, index was larger at the 
sampling site scale than at the community scale. 
Hence, the replacement of species between 
neighboring sampling sites is quite large in this 
short longitudinal range of the studied streams. 
In spite of this, the species richness, S, diversity 
parameter was statistically related between 
biological groups; e.g., fish species, FS, was 
negatively related to riparian tree species, PCS, 
(FS = 12.63 - 0.62PCS; 1?= 0.30; P < F = 0.001) 
as well as to benthic insect species, BIS 
(BI5.= 17,89. = 4,66 * La(EFS) 1*= 0,23; 
P < F = 0.024). Then, BIS is also negatively 
related to PCS (BIS = 17.83 - 4.66 * In (12.63 - 
0.062PCS). That is, the riparian tree community 
directly controls the diversity of the fish and 
benthic insect communities. Tree communities 
with a few typical riparian species colonizing 
upper river segments support the largest fish 
and benthic insect diversity and the strength of 
this relationship may fade away as lower river 
segments may have quite a small number of 
rare species. Lower river segments are the typi- 
cal habitat of Tilapia spp. and the relationship 
between riparian tree species and fish may shift 
during unseasonal wet episodes. 


Conclusions 


This research points to the finding of two 
distinctive riparian communities defined by 
the tree and fish communities in upper (RR, 
HCP, FR, ChP) and lower (MSJ, FSJ, PSJ, TR) 
river segments. The former cluster is de- 
scribed by pristine, high river discharge, cold 
river flows, and high water quality. The most 
distinctive tree and fish species are Taxodium 
mucronatum-Juglans spp.-Populus tremuloides; 
Notropis amabilis-Etheostoma grahami-Ictalurus 
lupus-Moxostoma congestum, respectively. Lower 
river segments are characterized by deterio- 
rated river flow quality, small discharge, and 
high flow temperatures. Tree and fish species 
assemblages colonizing these river segments 
are Cephalanthus occidentalis-Bacharis glutinosa- 
Prosopis laevigata-Fraxinus spp.-Celtis laevigata; 
Lepomis megalotis-Lepomis macrochirus-Gambusia 
affinis-Poecilia mexicana; respectively. Tree and 
fish communities show a strong within-group 
but a lack of between-group connectivity in up- 
per and lower sampling sites. The benthic insect 
community is classified into two clusters as well 
but they are more difuse indicating insects are 
better adapted to contrasting features of upper 
and lower river segments. It is concluded that 
riparian communites of Mexico's northern San 
Juan River are well adapted to variable river 
flow conditions controlled by hydro-climatic 
variability of several time scales. 
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